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We investigate the localization properties of the single particle spectrum of a one-dimensional 
speckle potential. We consider both the repulsive and the attractive cases. The system is controlled 
by a couple of parameters: (i) a dimensionless variable, determined by the mass of the particle, 
the correlation length and the average intensity of the field; (ii) the size of the box containing 
the particle. Depending on the value of these parameters and the considered energy level, the 
eigenstates exhibit different regimes of localization. We compared our numerical results to the 
present theoretical framework, showing the key role played by the finite size of the system and the 
self-correlation of the potential when the observed behaviour turns out to be different from the 
predicted one. 



I. INTRODUCTION 

Random potentials created by laser speckles are a pop- 
ular tool employed in experiments with ultracold atoms 
for investigating the behavior of disordered systems [1]. 
In particular, one-dimensional (ID) speckle potentials 
have been the object of an intensive study in recent 
years. Many interesting features have been addressed 
from both the experimental and the theoretical sides, in- 
cluding classical localization and fragmentation effects, 
frequency shifts and damping of collective excitations, 
inhibition of transport properties, Anderson localization 
and related phenomena [2HTq1. 

Despite this intense research activity, the properties 
of the single particle spectrum of the speckle potential 
have been addressed only partially, even in the ID case. 
For the uniform case, the low energy behavior has been 
first considered in ^ , where it has been shown that den- 
sity of states of repulsive speckles is characterized by a 
Lifshitz tail fTTJ [12] . This has been discussed more thor- 
oughly in [J^, for both attractive and repulsive speckles. 
There it has also been shown that, as for the charac- 
teristic extension and the localisation properties of the 
eigenstates, three different regimes of speckle intensities 
s can been identified: semiclassical {s ^ 1), interme- 
diate, and quantum {\s\ < 1). In the latter case, the 
presence of an effective mobility edge - separating local- 
ized and extended eigenstates - has been discussed in [9] , 
and its consequences have been observed in the Anderson 
localization experiment [10]. In addition, the existence of 
both extended and localized states in the presence of an 
inhomogeneous confinement has been discussed in [14 . 

What is lacking so far is a detailed numerical analysis 
of how the density profile of the eigenstates is related to 
the main properties of the potential (i.e.: spatial exten- 
sion, intensity, blue or red detuning), depending on the 
chosen energy level. 

In this article we consider a ID speckle potential de- 
limited by infinite walls - a sort of disorder in a box - 
and discuss the localisation properties of its eigenstates 
as a function of the speckle intensity s and the size L 
of the system. The box-like barriers can be realised ex- 
perimentally by focusing two laser beams perpendicular 



to the speckle direction, as discussed in [15 , allowing for 
the investigation of finite size effects in a textbook case. 

We studied both the ground state and the excited 
states, having the main target of characterizing the local- 
isation properties. In the two cases we used specific ob- 
servables, in order to explore the behaviour of the eigen- 
states while varying the parameters s (amplitude of the 
potential), L (size of the system) and the energy level 
E. When possible, we extrapolated the behaviour of the 
system in the limit L ^ oo. 

We observed a stochastic translation of the center of 
mass of some eigenstates while varying 5, depending on 
the single realisation of the speckle potential. Compar- 
ing a numerical measure of the probability of observing 
this phenomenon and a theoretical framework developed 
in this work, we separated the contribute linked to the 
boundary effects from the one depending only from the 
shape of the disordered potential. We expressed the lat- 
ter as a function of the partecipation ratio. 

According to the traditional literature, when L ^ oo 
a threshold energy level Ec{s) can be defined, that sep- 
arates localised from extended states (corresponding to 
an effective mobility edge). Moreover, this theoretical 
framework gives a closed form expression for the locali- 
sation length iioc as a function of ^ < Ec and s and a 
validity range for it. We performed closed form and nu- 
merical analysis of the validity range of the theory and re- 
stricted the interval of (5, E) where the predictions hold. 
Moreover, we described the behaviour of the system out 
of this new range in a closed form that fits our simula- 
tions. 

Starting from the traditional tools of multifractal anal- 
ysis, we defined a functional able to distinguish the lo- 
calised eigenstates from extended ones and from the ones 
belonging to the crossover region of the spectrum. We 
verified that this observable improves the performance 
of the traditional ones and we used it in order to char- 
acterize the three regions of the spectrum while varying 
s. 

The article is organized as follows. In Sec.|ll]we define 
both the considered system and the main tools that we 
used to characterize its eigenstates. Then in Sec. |III| we 
discuss the properties of the ground state and the dipen- 
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dence of its appearence from s and L. The density profile 
of the excited states is considered in Sec. |lVl Our results 
are summarized in Sec. [Vl 



II. MODEL AND METHODS 

Let us consider a single particle in a ID speckle po- 
tential Vs{x) = Vov{x/(^), with intensity Vb = (Vs) and 
autocorrelation length ^ [SJ [16]. Vq can be both pos- 
itive or negative, the potential resulting in a series of 
barriers or wells, that in the following will be referred as 
repulsive and attractive^ respectively. The autocorrela- 
tion ^ represents a natural length scale for the system, 
and = fi^/2m^^ the corresponding energy scale. The 
hamiltonian of the system can be written in a dimension- 
less form as 



H 



■ sv{x), 



(1) 



where the rescaled intensity s = 2m^^Fo/^^ represents 
the only parameter of the system. In case of a finite size 
system also the size L becomes relevant; its effects will 
be discussed hereinafter. 

The speckle pattern can be generated numerically as 
discussed in [5 (and references therein). Then, the 
spectrum can be obtained by solving the stationary 
Schrodinger equation HiIj = EtJ;^ by mapping on a 
grid with vanishing boundary conditions, corresponding 
to infinite walls. 

For a fixed system size L, the localisation behaviour 
of the eigenstates in the different intensity regimes of 
speckle potential have been outlined in [13]. The char- 
acteristic extension and the localisation properties of the 
eigenstates depends considerably on 5, and one can iden- 
tify three different regimes, ranging from the semiclas- 
sical limit, for |5| 1, down to the quantum regime, 
when 1^1 < 1 (with an intermediate regime in-between). 
In the case of attractive speckles, and for shallow po- 
tentials (quantum regime), the eigenstates extend over 
several potential wells. The lowest lying eigenstates lo- 
calise around deep wells but the ground state is not nec- 
essarily localised in the deepest well as the width of the 
well plays also a crucial role. When |s| is increased, the 
eigenstate width shrinks and their position may change. 
Eventually, for large enough \s\ (semiclassical regime), 
the eigenstates tend to stack up as bound states inside 
isolated wells, with the ground state being the lowest 
eigenstate of the deepest well. A similar scenario holds 
for repulsive speckles, in this case the localisation taking 
place in between barriers. 

Here we discuss it more thoroughly, characterizing the 
spectrum as a function of s, and discussing its depen- 
dence on L. In particular we consider a system length 
typically in the range L/^ G [100, 1000] (up to L = 7000 
in certain cases), and average over 50 to 500 speckle re- 
alisations (see Appendix [C] for further details). 



In order to characterise the localisation properties of 
the eigenfunctions we use the following quantities. 



1. The participation ratio 

PR [i/j] 



L 



J^dx\tlj{x)\ 



(2) 



that measures the relative extent occupied by the 
state ?/^; it decreases as gets more and more lo- 
calised. 

2. The localisation length iioc, that characterises the 
exponential decay of the tails of Anderson lo- 
calised states, |?/^(x)p oc ex.p{—2{x — xq) / iioc} ^ with 
xq = argmax | |?/^(x)|^| being the localization cen- 
ter. See Appendix [C] for further details on this 
quantity and the next two. 

3. The average curvature on the tails of the wave func- 
tion, defined through the following functional 



'd^ ln|V^(x)|^ 



dx'^ 



dx 



(3) 



where Li is the extent of left and right tails (over 
the ground noise), rj vanishes for exponentially lo- 
calised states, and is positive definite otherwise. 

4. The localization volume Dioc defined as the length 
of the interval over which the exponential localiza- 
tion occurs. 



5. 



The functional cr^ [•] (defined in ^IVB), that we 
introduce in order to characterize the eigenstates 
of the system by relating their localised/extended 
behaviour with their grade of self-similarity. 



The first two are widely used in the literature, while the 
others are defined here in order to show specific features 
of this system that more commonly employed functionals 
are not effective in measuring. 

In the next section we start by considering the prop- 
erties of the ground state as a function of 5, discussing 
their dependence on the finite size of the system. Then, 
in the following section we will consider the behaviour of 
the full spectrum. 



HI. THE GROUND STATE 

For a fixed system size L, the localization behavior 
of the eigenstates in the different intensity regimes of 
speckle potential have been outlined in ^3j. Here we 
discuss it more thoroughly. We start by considering the 
density distribution of the ground state for a single real- 
ization of length L = 400, as a function of s, as shown in 
Fig. [ij As an example, we discuss the case of a red de- 
tuned speckle realization; a similar behavior is observed 
for the blue detuned case as well. 
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FIG. 1. Behaviour of the density |'0(x)|^ of the ground state 
for a single reaUzation of red-detuned speckle potential (L = 
400). We have considered the following intervals of |s|, respec- 
tively from the top to the bottom: 2 • 10~^ < |s| < 2 • 10~^ 
(top); 4 < |s| < 40 (center); 10^ < |s| < 10^ (bottom). We 
plotted the considered values of |s| in increasing order from 
red to blue line. 



According to our simulations, some of the behaviours 
observed in Fig. [ijhave a general validity: 

{i) For too small values of |5| the localization behavior 
of the ground state is affected by the boundaries. In this 
case the tails decay exponentially, but the finite size of 
the system limits the localization length iioc to values of 
the order of the system size L. 

(ii) In the opposite limit, for large values of \s\ the 
ground state is bounded in a single well (red speckles) 
or between two adiacent barriers (blue speckles). In this 
case the tails ceases to decay exponentially. 

{iii) For intermediate values of |5| the behavior of the 
tails continuously connects the two extreme regimes. 

(iv) For some speckle realization, the localization cen- 
ter of the ground state may change as a function of s, see 
e.g. Fig. [2] This effect will be discussed in §111 B 

Two of the quantities defined in ^II]- 77(5) and Dioc{s) 
- allow for measuring the interval of |5| where there is 
a boundary effect (dark grey area), as shown in Fig. 
[3] for the same speckle realization used for Fig. [l] In 
fact Dioc{s) remains constant and r]{s) decreases until 
\s\ '> — 10~^, when the tails detach from the boundaries 
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FIG. 2. An example of the change of localization position for 
increasing \s\. In the example we considered L = 1000 and 
s ^ 0.1. The two ground states are obtained with a shift 
As 10"^ 

of the system. The decrease of 77(5) is linked to the lo- 
cal curvature of the tails when they are too near to the 
boundaries and, for even smaller \s\ values, to the fact 
that the ground state cannot present a completely expo- 
nential localization. When the boundary effects are over 
(Fig. [sj light grey area), there is an interval of |5| where 
r]{s) 0: di "pure" Anderson localization occours. By 
increasing \s\ we observe that the average curvature of 
the ground state begins to grow up, due to the increas- 
ing importance of the deep fluctuation of the potential 
where the speckle is localized. We have not observed any 
clear dependence of the amplitude of the region where a 
"pure" Anderson localization occours from s and L, since 
there is a very strong dependence of r]{s) from the shape 
of the single speckle potential, until it starts to constantly 
grow up. In fact, the most of the considered speckle real- 
izations show a noisy behavior of 77(5) for small \s\ values, 
due to the presence of big local fluctuations of the poten- 
tial and to random changes of the localization position. 

However, we observed on average an abrupt rise ofr]{s) 
for \s\ ranging from 1 to 3 (considering 100 < L < 500). 

A. Localization type 

We discuss here how the localization properties of the 
states - either exponentially localized by the full random 
profile of the potential or bound states of single fluctua- 
tions (in a single well or between two barriers) - depend 
on s and L. 

In the paragraph above we discussed the localization 
process of the ground state density |?/^(x)p at increasing 
1^1 . Fig. [4] describes the same process averaged over 50 
realizations of the speckle potential. We can observe that 
there is no transition from extended to localized states 
for both the attractive and the repulsive cases. In spite 
of this, both PR{s) and iLOc{s) exhibit a "knee" at in- 
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FIG. 3. ri{s) (blue dashed curve) and Dioc{s) (red solid curve) 
for the ground state of the speckle realization considered in 

fig- a 



creasing \s\. This knee separates the region where the 
localizing tails are influenced by the whole extension of 
the system (grey area) from the one where the leading 
contribute to the shape of |'^(x)p is given by the lo- 
cal configuration of the potential around the localization 
point (white areas). 

Thought the knees are not exactly coincident for 
PR{s^L) and £loc{s, L)^ their positions on s axis are 
comparable. The grey area extends conventionally from 
s —100 to s 300: in fact the localization process is 
faster for 5 < 0, according with the intuitive fact that a 
deeper well is more effective than two higher next barrier, 
in containing the localized ground state at increasing \s\. 
Moreover, the dipendence of PR{s^ L) and iLOc{s^ L) on 
L is neglectable (though present). As we will discuss in 
^ IV A[ we will see that this holds only for the ground state 
case, since L plays a more relevant role for the degree of 
localization of the eigenstates while increasing E. 

The curvature r]{s,L) is not less sensible to s varia- 
tions in the white areas than in the grey one. This is 
not surprising: even if an eigenfunction is so localized 
to be bounded in a single well, it goes on increasing its 
curvature as the depth of the well dilates. Moreover, 
r]{s^L) exhibits a non neglectable dipendence from L at 
1^1 ^ 1 values. This dipendence from the extension of 
the system is related to a local effect: in fact, a greater L 
implies a bigger number of wells (couple of barriers) and, 
consequently, a deeper global minimum (higher couple of 
next maxima) on average, where the ground state can be 
bounded more effectively and so exhibit a greater r]{s,L) 
value. The role of L is more relevant for 5 < 0, consis- 
tently with the greater binding effectiveness of the wells 
than the one of the barries, as stated above. 

In Fig. |5] we can observe the behaviour of Dioc{s^L) 
averaged over 50 realizations of the speckle potential. As 
stated before, Dioc{\s\ <C 1,L) ^ L. However, when 
1^1 > 1 the dipendence from L disappear for both s ^ 0. 
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FIG. 4. The functionals rj (top panel), PR (central panel) 
and Iloc (bottom panel) for the ground state as a function 
of s, increasing the size L of the system from L = 100 (red 
line) to L = 500 (purple Une) and L = 1000 (blue line). For 
each choice of s, L considered, the plotted value is averaged 
over 50 realizations of the speckle potential. 



B. Variation of the localization position 

In the following we analyze the phenomenon observed 
in figure |2j The change of localization position (CLP 
in the following) originates from the possibility that two 
next eigenstates exchange their position in the energy 
spectrum. This fact depends on the shape of the single 
realization of the potential and we found speckles without 
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cipation ratio PR{s, L) in a closed form. In fact we have 



Pcip L) = 



PR{s,L) 



27T 



(4) 



This result is obtained in Appendix |B] Our tests show 
that equation Q holds with a good degree of approx- 
imation. In figure [7] we plotted the crossover between 
two regions: the one where the boundary effects prevail 
[s 0) and equation Q is not verified and the other 
(large enough \s\) where the numerical simulation oipcip 
is well described by equation Q. 



FIG. 5. The functional Dioc for the ground state as a func- 
tion of s, increasing the size L of the system from L — 250 
(orange line) to L = 1000 (blue line). For each choice of s, L 
considered, the plotted value of Dioc(s, L) is averaged over 50 
realizations of the speckle potential. 



CLP events on all the tested interval of 5, as well as 
speckles which shows this behaviour one or more times. 

Given a set of speckle realizations and considering for 
each potential only the first CLP event observed for the 
ground state — where possible — while increasing \s\^ 
we build the cumulative density function Pcip{s^L). The 
result is shown in figure [6j 




FIG. 6. Probability distribution of observing one or more 
change of localization point event for lesser than or equal 
to a given \s\. The function was obtained simulating the be- 
haviour of the ground state over 500 speckle realizations as a 
function of |s|, for both red detuned (top) and blue detuned 
(bottom) speckle potential. 
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FIG. 7. Probability density function for the occurrence of 
a CLP event. Red circles are computed from the numerical 
results plotted in figure [6] for L = 1000; the blue solid line is 
obtained from equation l4| . 



IV. EXCITED STATES 

Let us now consider the behavior of the excited states 
as a function of their energy E. We expect the lower part 
of the spectrum to be affected mainly by the presence of 
the "disordered" potential, whereas at high energy the 
finite-size effect (the effect of the boundary conditions) 
should be predominant. As an example of the behaviours 
discussed in thi paragraph, in Fig. [s] and |9] (top panels) 
we show a density plot of the eigenstates, for a single re- 
alization of the speckle potential, for both attractive and 
repulsive cases ( the fu nctional and the lower panels 
are discussed in pV B ) . We observe the presence of three 
regions: localized states for low energy, extended states 
for high energy, separated by a wide intermediate region. 
In the following we analyze the localization properties of 
the eigenstates as a function of their energy by consid- 
ering their multifractal behavior and localization length. 



With the exception of a small interval of |s| :^ 0, where 
the boundary effects play a major role in the energy level 
of the ground state, the probability density function pdp 
of observing a CLP event can be linked with the parte- 



A. Localization length 

In |9] is shown that, under the Born approximation, 
there is a high-momentum cutoff kc (that represents an 
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effective ME)^ such that ah the matter waves with mo- 
mentum k > kc are extended on every length scale (i.e. 
ilhcik > kc) = 0). 

In our case (see also [3 [T3^) the expression for £loc 
given in [9_ reads (see Appendix [A]) 



7.04 



E- 



with 



-oo 



Er 



E <Ec 
E>Ec 



(O.SStt)^ 



and under the condition 

\s\ < sc = 2(0.887r)2 . 



(5) 



(6) 



(7) 



The Mobility Edge is not directly observable using nu- 
merical tecniques. In fact, an infinitely long tail requires 
an infinitely extended system in order to be free from 
boundary effects. However, from equation ^ we know 
that the following should be true 



lim £iocis,L,Ec (s)) 

L^+oo 



+00 



(8) 



for all s that satisfy equation So, a numerical marker 
of the presence of a Mobility Edge is that £ioc{s^ Eq) 
is increasing, as a function of L. In order to investigate 
the presence of this marker, we simulated the eigenstates 
of the system around Ec{s)^ using different values of s 
and the largest values of L we can numerically access at 
the moment. As shown in figure [lOj we can observe the 
following outcomes 
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FIG. 8. (top) Example (s = 1) of density plot for the spec- 
trum of a red detuned speckle potential, (center) Compari- 
son between the eigenstates ipE of a, single realization of the 
speckle potential and [^Pe]- This will be discussed in ^IVB 
(bottom) Density profiles of three eigenstates are plotted in a 
convenient logscale, in order to show the differences amongst 
the localized, delocalizing and delocalized wave functions. 



• There is a range of values of s that verifies equation 
([8| , suggesting the presence of an effective Mobility 
Edge in the limit L ^ oo. 

• The boundaries of this range are not the same for 
the red— detuned and the blue— detuned potentials. 
Equation ^ becomes 



' < 5 < 5 



blue 
C 



where -2 < s^J"^ < -1 and 0.1 < 



< 1. 



(9) 



and 5^^^ have not to be intended as critical val- 
ues of a transition between two regimes (Mobility 
Edge/other), because the condition to be satisfied 
is an approximation. Considering this, we can ex- 
pect that the passage from the Mobility Edge oc- 
cours with a crossover. 

• For both the attractive and the repulsive potential 
and for L 1, iioc seems to be indipendent from 

red r>blue1 



L, considering L > 5000 and s ^ [s^^, s 



c 



From figure 10 we can observe that iLOc{s = ±0.1, L > 
4000, E = EcJ— L. This means that the eigenstates are 
extended for the tested values of L. However, the fact 
that the eigenstates go on being delocalized at £^ Ec 
for increasing values of L is compatible with equation Q . 

For lesser values of L than the range tested above, we 
could measure iLOc{E) starting from the ground state 
to £^ = Ec and beyond. Until Iloc is comparable to 
L, we can test the prediction of equation ([5| against the 
simulation being free from boundary effects. 

As shown in figure [TT] the numerical measure of 
^LOc{s^ E) fits the the trend predicted in equation ([5| 
only for s = {±0.1, —1} (until the extension of the tails 
reach the boundaries of the system). The same values 
of s satisfy equation ([8|, confirming the range of validity 
estimated in equationT9|. We have seen above that Iloc 
becomes indipendent from the size of the system for large 
enough values of L, for blue-detuned potentials out of the 
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FIG. 9. This figure is analogous to Fig. |8]for a blue detuned 
speckle potential. In this case we have considered s = 1. 
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FIG. 11. iLOc{E) for -10 < s < 10 and L = 800. The 
figure shows a comparison between the simulations (orange 
dots, averaged over 50 realizations of the potential) and the 
theoretical results for L = +oo (continuous line). 



FIG. 10. Localization length iioc {s, L, E) measured for L G 
[4000,7000] and E G [Ec{s) - 0.25, Ec{s) + 0.25]. For each 
value of L and s, iioc is averaged over the range of E cosidered 
for 50 realizations. 



range stated in equation ([9|. In the following, we study 
the approach to -^iLoc{s, E) = for increasing L. 

The figure 12 confirms the non comparability between 
the prediction of equation ([5| and the numerically ob- 
served behaviour of iLOc{s = l,L,E), considering L G 
[200,900]. In order to describe this behaviour, we fol- 



lowed the steps shown in figures 13 and 14 



1. We looked for the function of E which fits the simu- 
lations better, keeping s and L constant. As shown 
in figure [13) we have: 
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hoc{E)U^L = exp{(/i \n{E) + 52} 



(10) 
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FIG. 12. Iloc{E) for s = 1 and L = 200 (orange dots), 
L = 900 (blue dots). The figure shows a comparison between 
the numerical simulations (dots, averaged over 50 realizations 
of the potential) and the theoretical prediction (continuous 
line). The dashed line marks Ec- 



2. We looked for the functions of L which fit gi and 
g2 better. As shown in figure [Ml we have: 



5i(L) = aiL/(/3i+L) 

52(i) = a2/(/32 + i)+72 



(11) 



which imply the asympotic indipendence of Iloc 
from L in the limit L ^ oo. 



Generalizing the results contained in equations ( 10 ) and 
(11) for all s > s^^^, we can suppose that 



W(S,i,i?).>.«"e (Xi?''(^) (12) 

for L ^ 1 (expressed in ^ units). Testing several s 



values, we found that the function h{s) in equation (12) 
which best fits the results of our numerical simulations is 
h{s) = /iil^l^^ xhe same approach described above can 
be applied to the red-detuned case, leading to the same 
results. Finally we have: 



^LOC 



(s,L,^) = /io^^^''''' 



Eq. 



5000, 



^ ^gi^d gblue I 



(13) 



and 



(13| holds for L > 

E > max(0,£^^s). Table IV A reports an estimate of the 
parameters ho, hi and /i2- 
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FIG. 13. iLOc{E) for s = 1 and L between 200 (orange) and 
900 (blue) . The top figure shows the numerical dipendence on 
L while increasing E. Iloc does not show any sensibility to 
the ME (dashed line) . The numerical results are averaged over 
several realizations of the potential. The bottom figure shows 
a fit of the same output, using the function Iloc{E^ L) — 
exp{^i(L) hi{E)+g2{L)} (chosen on the basis of the numerical 
evidence) . The average value of the of the fit is 0.985 across 
the fitted functions. 
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FIG. 14. The figure shows the fitted values of gi and g2 
(dots) across the tested values of L, plotted together their 
respective fit errors (bars), considering s = 1 as an example. 
The numerical results are fitted using the functions described 
in equations (10) and (11). 



Self similarity of eigenstates as a marker of the 
delocalization 



For finite size systems with a J— correlated disordered 
potential, the presence of a transition from localized 
to extended states has been found both for the three- 



sign{s) 


ho 


hi 


h2 


+ 


0.590 


2.927 


-0.901 




5.048 


2.209 


-0.706 



TABLE I. Best estimate of the parameters used in Eq. |13| 



dimensional [17] and the one-dimensional case [18] (the 
so-called Anderson transition). For these systems, the 
eigenstates around the critical energy threshold show a 
kind of scale invariance known as mult ifr act ality. 

Generally speaking, a scale invariance property of a 
given function f{x) can be numerically verified partition- 
ing it into clusters of increasing size. The necessary con- 
dition that f{x) must satisfy for being considered as a 
fractal (or a multifractal) is the independence of a cer- 
tain quantity uj (g. A) - defined later on - from the cho- 
sen size A of the clusters. In the following we specialize 
f{x) = |?/^(x)p. Let us introduce the measures 



(A)= / 

JkA 



(fc+l)A 



dxf{x) k = 0,...,kmax{^) (14) 



where the interval [0, {kmax + 1) A = L] is the spatial 
domain of the system for every kmax ^ N. For each q eZ 
and for each acceptable value of A, we can compute the 
q—th moment Pq{A.): 



For a self-similar function, it holds that 
P,(A) (X A'^(«). 



(15) 



(16) 



When the relation (16) is verified, the shape of uj{q) de- 
fines the type of self-similarity of the function (if uj{q) 
is linear with a single slope, f{x) is fractal; if uj{q) is 
linear with two different slopes for positive and negative 
(7, f{x) is multifractal; otherwise f{x) is not a fractal in 
any sense). For every f{x) (even for the non self similar 
ones), the following quantity is well defined: 



j{q,A) = lnP<^(A)/lnA. 



(17) 



Given that f{x) can be a (multi) fractal only if the nu- 
merical dipendence of cj(g,A) from A is almost absent, 
we can introduce a quantity to measure how much far is 
f{x) from being self-similar in some sense: 



[f]=^U{q.A)-{uj{q.A))l 



(18) 



The results obtained in [TF suggest that, if we can ignore 
the effect of the correlation in the speckle potential, our 
system should show an Anderson transition at a critical 
energy Ec^ with a multifractal behaviour of the eigen- 
states with eigenvalues near to Ec- This is in agreement 



with the prediction of a ME in ^IVA and it should lead 
to observe numerically 



lim 

E^Ec 







(19) 



We notice that both the predictions of a ME stated in 



^IVA and in this paragraph do not take into account the 
effect of the spatial correlation of the speckles at high 
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orders. In particular, the presence of a multifractal be- 
haviour has been observed in completely decorrelated ID 
system. 

Looking for a quantity able to distinguish among dif- 
ferent degree of localization, we compared the results ob- 
tained with [i/je] against the other main observables 
defined in qlL and the variance of the density H^e]- 
This functional turns out to be the only one able to react 
at the crossover between localized and extended states. 

The link between the possible density profile of the 
wafefunctions and the correspondent values assumed by 
[iI^e] is shown in Fig. [s] and [oj For both localized 



and extended states, [ipEl oscillates randomly around 
a constant value. In the intermediate delocalization re- 
gion of the spectrum, [i^e] decreases while E is grow- 



mg. In order to measure and E^^ (respectively the 
attachment point and the detachment point of the inter- 
mediate energy tranche where the delocalization process 
occours), we fitted Incr^ H^e] with a polygonal chain (see 



Appendix |C| for further details). 



The behavior of E, 



(a) 



and E'c' as a function of s is shown in Fig. 



15 



We 



have also checked that the extent of the intermediate re- 
gion is not affected by the system size L for small enough 
systems {L < 10^). 



- (L= 500) 




fully extended eigenstates 






— fully localized eigenstates 







FIG. 15. The amplitude of the intermediate region as a func- 
tion of s. E^^^ an d E^^ ^ are measured averaging the values 
obtained from eq. (CI) over 50 realizations. The dashed line 



plots the predicted threshold for the Mobility Edge. 



Fig. 15 shows the localizing behavior as a function of 



• For 5 > 0, the intermediate region is shifted to- 
wards higher energy levels and its amplitude in- 
creases, for increasing values of s. This fact can 
be related to the distribution of the maxima of the 
potential, that are higher and more dispersed for 
larger values of s. Conversely, for s < the inter- 
mediate region is far more stable in mean and am- 
plitude, according with a distribution of the max- 
ima less sensibles to the s variations and more con- 
centrate around its mean value. 

• The prediction for the effective ME Ec{s) in eq. (|6| 
is not in agreement with the position of the inter- 
mediate region [E^^\eP]. This is coherent with 



the results discussed in ^ IV A , where we verify that 
the finiteness of the system forbid the numerical 
observation of the predicted Mobility Edge. 



• Eq. (19) is not verified. This can be related to the 
presence of self-correlation in our potential (unlike 
the case studied in [18 ). Once again, nothing sug- 
gests the presence of an Anderson transition for the 
considered system. 

The last two points suggest that the finite size of the 
system and the correlations in the disorder play a crucial 
role in the properties of the spectrum generated by a 
speckle potential and they cannot be neglected at any 
order. 



V. SUMMARY 

In this work we focused on the density profile of the 
eigenfunctions associated to a particle in a box-plus- 
speckle-potential system. We analyzed how the parame- 
ters 5, L and E infiuence the density profile of the eigen- 
functions, giving special attention to the degree of local- 
ization. 

The ground state exhibits a very poor dependence from 
L. The passage from Anderson localized states to eigen- 
states bounded in a single well (or couple of barriers) is 
smooth at increasing \s\ but it is stil possible to partition 
the localization process into two distinct regimes. The lo- 
calization is faster (as a function of \s\) for the attractive 
case. 

We found that the lower energy eigenstates can ex- 
hibit an abrupt shift of the localization point (and so of 
the center of mass of the system) while tuning due 
to an occasional degeneracy with the next energy level. 
We linked the probability of observing this effect to the 
Partecipation Ratio of the eigenstate and to the finite 
extension of the system - the latter contribute being rel- 
evant only for |5| <C 1. 

The system does not exhibit any Mobility Edge nor 
multifractal behaviour while delocalizing at increasing E. 
The first fact is attributable to the finite extension of the 
system, while the second is linked to the self-correlation 
of the potential. We described in a simple closed form 
the asymptotical behaviour of the localization length as 
a function of s and for L 1 and for s out of the 
validity interval of the mobility edge framework. Finally, 
we found that the eigenstates of this system at increasing 
energies can be divided in three distinct regions according 
to their degree of localization. The quantity developed 
here in order to observe this partition could find an appli- 
cation also in other contexts as an indicator of the degree 
of self-similarity in data series. 
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Appendix A: Generalization of the ME to finite-size 
systems 

As stated in [9 , the Born approximation in this context 
assumes the inequality 



(Al) 



where ^ = O.SStt^. Using in Eq. (Al) the definition of 
s given in Section [llj we obtain the following equivalent 
condition: 

E^s^Cs^^^, (A2) 
0.2. Under this condi- 



where C = 2-4/3(0.887r)-2/3 
tion, we have 



-m 



(A3) 



where ^JVi2j(k) = ^/^TJ2{l-kl)Q(l-ki) is the autocor- 
relation function of the potential in case of an infinitely 
extended system. Substituting the definition of s again, 
Eq. ( A3 ) in our framework becomes 



1 



(A4) 



which leads to eq. ([5|. From eq. (A2) and (A4), it 
follows that 



s + Cs^^^ <^E <s^ (0.887r)2 



(A5) 



that leads to a stronger requirement than ( Al ) for allow- 
ing the good definition of Iloc^ as stated in eq. 



As shown before, low energy eigenstates of the discussed 
model are well described by 



(B2) 



where Dtaiis is the domain where the wavefunction ex- 
ponentially decays and g{x) is the ground noise out from 
Dtaiis- We want to explicit the value of Ei^ in order to 
solve the equation (Bl). We can state that 



E,, 



(B3) 



using the fact that the hamiltonian is separable. For the 

kmetic term, approximatmg g{x) {) ^ e ^loc , we 
have 



L g2 



^LOC 

1 



Jo 



^Loc dx 



(B4) 



this approximation holds for \s\ large enough that bound- 
ary effects are neglectable. For the potential term we 
have 

= r \M^)\\{x)dx ^Y.^lvj (B5) 
Jo ^ 

where we approximate the continuos system applyng an 
arbitrary discretization. Neglecting the self correlation of 
the potential, we can consider all the Vj as i.i.d. stochas- 
tic variables. Since Vj follows a negative exponential dis- 
tribution, we have 



(B6) 



A sum of indipendent exponentially distributed variables 
with different mean life values fits the generalized Erlang 
distribution. So we are able to determine its variance in 
a closed form. 



(B7) 



Appendix B: CLP effect 

In the following we are going to prove that the result 
in equation Q holds except for a small interval around 
5 = 0. Naively speaking, a change of localization point 
event (CLP event in the following) occours when two 
next localized eigenstates exchange their positions in the 
spectrum because of a variation of s. Without loss of gen- 
erality, let us consider the ground state ipo and the first 
excited state ipi , with their respective energy eigenvalues 
^0 and El. So we have 



PC LP {s) =p{Eo = Ei\s). 



(Bl) 



On average, it holds that cr^ [Vq] = cr^ [Vi] and /i [Vq] = 
/i [Vi]. Approximating the distributions of Vq and Vi with 
two normal distributions, we have 



£^0 — ^1 



ms.L)^_._^ (B8) 
27r 



Rep lacing equations (B3), (B4) and (B8) into equation 
(Bl ), we obtain 



PCLp{s, L) = p{v = 0|5, L) 
that leads to equation Q. 



PR{s,L) 



27T 



(B9) 
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Appendix C: Numerical methods and criteria 

In the fohowing section we summarize some methods 
that we used to obtain the results contained in Sec. IIIII 
andHVl 

Generation of the eigenstates. We use the follow- 
ing prescription: Hij = 1/Ax^ {2dij — — ^i-ij) + 
sI{jAx)Sij^ with vanishing boundary conditions, ipo = 
V^AT+i = 0. The number of grid points N is chosen so 
that every fluctuation of the potential is sampled with at 
least 8 points (i.e.: N > 8L/1). 

Identification and quality check of the exponential tails. 
We need to identify four x values (xi,...,X4 in as- 
cending order) in order to select the two tails prop- 
erly. The routine starts choosing the points X2 and X3, 
near the localization point xloc = argmax[|?/^ (x) \] but 
far enough to exclude the local fluctuation of the po- 
tential around which the eigenfunction is localized (i.e.: 
X2,3 = xloc^O- 

The terminal point of the right tail x^ is found going 
through the |'0(^)| values for increasing x > xs^ until 
< ipnoise OT X > L — The routiuc identifies Xi 
using a symmetrical rule to the one described above. The 
threshold t/Jnoise has the form 

where a is a user defined parameter. 
We checked the exponential shape of the tails system- 
atically, computing the i?^ for the fit used to measure 
^LOC- When R'^ value returned an ambiguous indication, 
we verified the shape of the tails by a direct observation. 

Check for the physical meaning of the ground noise. 
In §[TT]we introduced the localization volume Dlqc to 
better understand the local influence of the boundaries 



over the density profile. Dlqc is well defined only if 
the abrupt interruption of the exponential decay of the 
eigenfunctions out of the interval [xi,X4] is not a nu- 
merical effect - at least in the discussed region of s, L. 
We verified this fact performing the three tests described 
below over a sample of speckle realizations: {%) computa- 
tion of the same eigenstate using different computers; (n) 
computation of the same eigenstate using a re-sampled 
version of the speckle potential with a different N (pro- 
duced by a spline interpolation); {Hi) comparison of the 
stability of iIj{x) to an "adjusted version" of the same 
eigenstate il){x). The latter is obtained extrapolating the 
esponential decay outside [xi,X4] until the machine dou- 
ble epsilon (~ 5 • 10~^^^ in this case) or a boundary of the 
system is reached. This test was performed applying the 
following condition to check the considered eigenstates 



< 



where fi\v ^ R^] := v/\v\. Ideally it should hold that 

n o I-L^Ije = I-L^I^e/ but we prefered the criterion above 
to avoid the error over the estin jate cj f the eigenvalue E. 

Measure of E\ 
in Fig. ^ M and 



(a) 
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and E^^^ in \ 
are obtainec 



IVB 



The results shown 
Tittmg In [crl] {E) with 
a polygonaT chain. More in detail, we used the following 
function 

&l [^e] = exp{(ci^ + C2) e{E^c'^ - E)e{E - Eg^) 
+ {csE + C4) 6(^ - eI^^)G{E^c^ - E) 

+ C56(^-4'^)} 

from which we obtained an estimate of E^^ and E^\ 
respectively the attachment and the detachment point of 
the crossover region. 
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